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Abstract. In most configurations aimed at generating entangled photons based on 
spontaneous parametric down conversion (SPDC), the generated pairs of photons are 
required to be entangled in only one degree of freedom. Any distinguishing information 
coming from the other degrees of freedom that characterize the photon should be 
suppressed to avoid correlations with the degree of freedom of interest. However, this 
suppression is not always possible. Here, we show how the frequency information 
available affects the purity of the two-photon state in space, revealing a correlation 
between the frequency and the space degrees of freedom. This correlation should be 
taken into account to calculate the total amount of entanglement between the photons. 
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1. Introduction 

Pairs of photons generated via spontaneous parametric down conversion (SPDC) are 
widely used in quantum cryptography, to test violation of Bell's inequalities and in 
other quantum information tasks, including dense coding and teleportation pQ. The 
popularity of this source of paired photons is strongly related to the relative simplicity of 
its experimental realization, and to the variety of quantum features that down converted 
photons can exhibit. It is well known that a pair of photons generated via SPDC can be 
entangled in different degrees of freedom, for example in polarization [21 E], in frequency 
[HE] and in the equivalent degrees of freedom: orbital angular momentum j6j[7j, space 
and transverse momentum [8j [9] . 

Usually, when considering entanglement between the two generated photons, only 
one of the degrees of freedom in which the pair may be entangled is considered and 
any influence of the others is neglected. This is often justified by the assumption of 
a filtering process in which the use of very narrow filters in one degree of freedom 
allows to neglect it [10]. Only one degree of freedom is considered, for example, in the 
experimental demonstration of teleportation [H] and in some proposals to measure the 
spatial entanglement between paired photons [8]. 

Pairs of photons generated with specific characteristics in various degrees of freedom 
have been also considered. For example, in hyperentanglement [T2l [T3] where the 
photons are entangled in various degrees of freedom that are assumed to be independent 
and, in configurations where one degree of freedom modifies the quantum state of the 
photons in other degrees of freedom: to distillate position entanglement using the 
polarization [H] or to control the joint frequency distribution of the generated pair 
by modifying the spatial properties of the pump beam [5]. 

All the configurations mentioned before assume a specific relationship between 
frequency, polarization and space in the two-photon state. In this paper, we will 
concentrate on the frequency and spatial properties of a pure two-photon state. We 
will show that the purity of the spatial two-photon state generated by SPDC depends 
on the frequency information available. This imply a correlation between space and 
frequency that cannot be ignored in any real filtering process, where the finite bandwidth 
of the filters affects the separability of the frequency and the space degrees of freedom. 
This is analogous to the effect of the spatial [15] and frequency [16] correlations of the 
paired photons on the degree of polarization entanglement of the generated photons. 
Additionally, in this letter we will show how the spatiotemporal correlations affect the 
entanglement between the photon. We will calculate how the purity of the signal photon 
in space and frequency, and therefore the amount of entanglement between the photons, 
depends on the geometry of the SPDC configuration, the characteristics of the pump 
beam and the filtering process. 

This paper is divided in three sections. In section [2] we discuss the main 
characteristics of the two-photon quantum state generated in spontaneous parametric 
down conversion considering the space and frequency degrees of freedom. In section El 
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the spatial state of the two-photon is described. It is shown how the spatial state may 
become mixed as a consequence of tracing out the frequency variables, which implies 
correlation between those degrees of freedom. Finally, in section HJ the purity of the 
signal photon state is analyzed taking into account the space and the frequency degrees 
of freedom. 



2. Quantum State of pairs of photons generated via spontaneous 
parametric down conversion 

2.1. General case 

Spontaneous parametric down conversion (SPDC) is an optical process in which a 
nonlinear crystal of length L is illuminated by a laser pump beam propagating in the 
z direction. Due to the interaction of a pump photon with the crystal there is a small 
probability to generate a pair of photons: signal and idler. In general, the direction of 
propagation of the pump beam is modified by a small angle po due to the presence of 
Poynting vector walk-off inside the crystal. All photons that interacts in the nonlinear 
process are characterized by the transverse momentum q n = (q%, q%), and the frecuency 
uj n = u® + Q n , where is the central frequency and fl n is the frequency deviation. The 
index n labels the signal (n = s), idler (n = i) and pump (n = p). 

For the sake of simplicity, in what follows we will assume that in the frequency and 
space degrees of freedom, the two-photon state is pure [17] , therefore the density matrix 
reads 



dq s dQ s dq^idQ i dq' s dQ,' s dc^ i dQ,' i 

x$(q s ,^,q 4 ,a)$*(q;,^,q;,^; 

x |q s ,fi s )| qi ,a)(qlX|(q-,^-|. 



(1) 



The spatiotemporal mode function $(q s , fl s , q^, fli) in a noncollinear configuration is 
given by 



$(q s , n s , q 4 , Qi) (X E p (q x s + gf , A )B p (Q a + Q,) (2) 

X C spatial ( q^ ) C spatial ( qi ) 

X J~ frequencyiS^s)^ f requency iS^i) 

/A k L s 

x sine 



where B p (ojp + Q p ) and E p (c[ p ) are the frequency and transverse momentum distribution 
of the pump beam at the center of the nonlinear crystal, respectively. The effect of 
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the unavoidable spatial filtering produced by the specific optical detection system is 
described by the spatial collection function C spat i a i((in)- The presence of frequency filters 
in the paths of signal and idler is indicated by frequency filter function J r f re quency(Qn)- 
The Phase matching conditions appears in equation ([2]) through the "Delta" factors, 
defined by 



A = q y s cos ifs + qf cos ifi + k s sin ip s - ki sin ifi (3) 
Afc = k p — k s cos ip s — ki cos <fi — q y s sin ip s + q\ sin ifi 
+ (q* + ) tan p cos a + A tan p sin a 



where k n = 
and fi are t 



[uj° n n n /cY 



1 1/2 



|q« 



is the longitudinal wave vector inside the crystal. ip s 



re propagation directions of the generated photons inside the crystal with 



respect to the pump direction z and a is the azimuthal angle [TOj [18] . 



2.2. Gaussian approximations 

In order to find analytical expressions for the purity of the spatial two-photon state and 
the purity of the signal state in frequency and space, we make some assumptions and 
approximations following refs |19[ [20] . 

We consider a gaussian shape for the spatial and frequency distributions of the 
pump beam: 



Epilpi Qp) = exp 



w 

-((? 

4 vy P 



x2 



and 



B p (Q s + Qi) = exp 



T 2 
4 p 



(4) 



(5) 



The pump beam waist is given by w p and the pulse temporal duration is T . 

Both, the spatial collection and the frequency filter functions are described 
by Gaussians, so that C spatia i (q n ) oc exp [-w 2 |q n | 2 /2] and T frequency (Vt n ) oc 
exp [— Q 2 l /(2B 2 )} where w n is the spatial collection mode width and B n is the frequency 
bandwidth. Notice that, in order to have more convenient units, the filters in momentum 
are defined in a different way than the filters in frequency. While B s i — > implies a 
single frequency collection, the condition for a single q vector collection is w s> i — > oo. In 
what follows we consider B s = Bi and w s = w^. 

Under these considerations and by making the approximation sinc(AfcL/2) w 
exp [— /9 2 A 2 L 2 /4] with f3 = 0.455, the mode function of the generated pair can be 
written as 



$(q s ,fi s ,qi,fii) = iVexp 



(6) 
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where N is a normalization constant, A is a positive-definite real 6x6 matrix containing 
all the parameters that describes the SPDC process and x l is the transpose of the vector 
x given by x t = ( q*, q y s , Q s , qf, qf, 

The physical system described by equation (JBJ) can be considered as composed by 
two subsystems. These subsystems can either be the two degrees of freedom, figure Q] 
(a) or the two photons, figure [I] (b). In both cases, the correlations that may exist 
between the subsystems imply that the system is not in a separable state, i.e., when 
considering the two degrees of freedom $(q s , fi s , q^, fij) ^ < & q (q s , q;)^^? ^i) and 
when considering the two photons $(q s , fl s , q^, fij) 7^ $ s (q s , f2 s )$j(qj, fij). 



(a) 



(b) 
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Figure 1. The two-photon state in space and frequency can be considered as composed 
by two subsystems. In (a) those subsystems are the degrees of freedom. The frequency, 
in gray, is traced out in the calculations of Sec. [3] In (b) the subsystems considered 
are the generated photons. In gray we show the idler photon that is traced out in the 
calculations of Sec. 0J 



3. Correlations between space and frequency degrees of freedom 

Let us consider the two degrees of freedom as subsystems like in figure [T](a). The state 
of one subsystem can be calculated tracing out the parameters that describe the other 
subsystem. Tracing out the frequency, the reduced density matrix for the resulting 
spatial two-photon state is given by 



P q = Tr n (p) 

dq s dO s dqj jdq' s (iq^ 
$(q s , Q s , q h fii)$*(q^ fi„ q-, ^)|q s , qi)(q' s , q-| 



(7) 



and the corresponding purity Tr(p^) of the spatial two-photon state is 
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x $(q s , n s , q i; ^)$*(q' s , ft 8J q-, ^ 

x$( q ;,fi;, q :,^)$*(q s ,^,q,,fi- 



(8) 



We can solve this integral using the exponential character of the mode function given 
by equation so that the purity becomes 

m , 9N det(2A) . , 

Tr(p 2 a ) = . K ' 9 

where B is a positive-definite real 12 x 12 matrix given by 



N exp 



1 



X l BX 



$(q s ,fi s ,qi,fii 



(10) 



X* = ^ x*, x'* J is a vector resulting from the concatenation of and x' 1 = 

(gf, gjf, Q' s , qf, qf , 

Since the composed quantum system described by equation (J2j) is in a pure state 
(Tr(p 2 ) = 1), the purity of the spatial two-photon state given by equation [91 is different 
from 1 only if the degrees of freedom are correlated [21]. This correlation is stronger 
when the partial purity gets closer to P . 

Notice that the purity of the spatial two-photon state is equal to the purity of the 
frequency two-photon state when the spatial degree of freedom is traced out. Therefore 
all the results derived for the spatial two-photon state can be extended to a frequency 
two-photon state. 

We plot the purity of the spatial state Tr(p^) as a function of the spatial filter 
width for different frequency filter bandwidths, figure [2j The half width at 1/e of the 
frequency filters in the wavelength variable is given by AA n , that is related to B s and 
Bi by B n = ircAX n / (A 2 Vhi 2). The SPDC process considered occurs in a L = 1 mm 
lithium iodate (LHO3) type I crystal, illuminated by a pump beam with A° = 405 nm 
and a beam waist w p = 400 /mi. The generated signal and idler photons are emitted 
with a wavelength X° s = A° = 810 nm at an angle (p = <p s = ipi = 10°. We consider 
negligible Poynting vector walk-off (po = 0) and in this case all the values of a are 
equivalent. 

In figure [2] we see that the purity of the spatial two photon state can be different 
from 1, which implies a correlation between space and frequency. The strength of the 
correlation increases as the purity decreases. From this figure we can see that as w s 
and Wj increase the spatial purity of the two-photon state gets closer to 1. This can be 
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Figure 2. Purity of the spatial state of the pair of generated photons for different 
values of the frequency filter bandwidth, as a function of the spatial collection mode 
width (w s = Wj). The pump beam waist is w p = AOOfim and the emission angle is 
<p = 10°. 



understood intuitively since an infinitely big w s (wj) implies the case of collecting only 
one q s ( q$) and therefore the two-photon state is separable in frequency and space. 

The separability and the lack of correlation between frequency and space, can also 
be seen in the "narrow band" limit where AX S = AAj — > nm: for any value of w s = w$ 
the purity is always 1. For a typical commercial interference filter, AX S = AAj ~ 1 
nm, the correlation between space and frequency is significant and consequently the 
purity decreases as can be observed in the figure. As the width of the frequency filters 
increases, the purity converge quickly enough to make the case of AA S = AA.j = 10 nm 
almost equivalent to the case without frequency filters AA S , AAj — > oo. 

In order to compare these theoretical results with previous works, in figure 3 we 
show the purity of the spatial state as a function of the collection mode, for three 
reported experiments. In figure [3] (a) we used the data of reference [5] . In their case, a 
type I SPDC process in a 1mm long LHO3 crystal is analyzed. As a pump beam their 
considered a diode laser with wavelength X p = 405nm and bandwidth AA P = 0.4nm. 
The pairs of photons generated propagate at (p = 17° and are detected by using 
monochromators with A\ Sti = 0.2nm. Two cases are depicted, the first one when 
the pump beam is focalize with a waist of w p = 30/im over the crystal and the second 
one where the pump beam waist is w p = 462/im. Their main result is the control of 
the frequency correlations by changing the spatial properties of the pump beam. As the 
pump beam influence the spatial shape of the generated photons, they report implicitly 
a correlation between the space and the frequency. For their spatial collection mode 
w s = 133, 48/iin it can be seen in the figure that according with our theoretical model 
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Figure 3. Spatial purity as function of the collection mode for three experimental 
cases. In (a) is showed the case of reference [5], in (b) the case described in reference 
[22] and in (c) the case reported in reference [23] , 

the correlation appears. 

Figure 3(b) shows Tr[p^] for the case described in reference [22]. A continues wave 
pump beam with X p = 405nm is used to generate photons with X si = 810nm in a 
1.5mm length BBO crystal. Interference filters of AA Si j = lOnm are used before the 
detection. The pump beam waist satisfies the condition w p » L and the generated 
photons propagates at tp — 0°. As it can be seen from the plot, the lack of correlation 
in this case is due to the collinear configuration. Under the conditions reported, but 
in a non collinear configuration will not be possible to neglect the space and frequency 
correlations. 

Finally in figure 3(c), Tr[/?q] is plotted in the case described in reference [23]. In 
their case a pump beam with X p = 351. lnm is focused with a waist of w p ~ 20fim on 
a 2mm long BBO crystal. Generated photons at X Sj i = 702nm, propagating at 4°, are 
collected after passing lOnm interference filters. In this noncollinear case, the lack of 
correlations are due the small pump beam waist. As can be seen in the figure, for less 
focused pump beams the correlations become appreciable. 

4. Spatiotemporal entanglement between signal and idler 

The correlations between frequency and space in the two photon states showed in the 
previous section, suggest that for a complete description of the entanglement between 
the photons, both degrees of freedom need to be taken into account. 
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In contrast with the previous section, here the physical system considered is 
composed of two photons, each one described with space and time variables, like is 
shown in figure [T] (b). In what follows, we will calculate the purity of the signal photon. 
Since the global state given by equation (|2J) is pure, the result of this section allows us 
to calculate the degree of spatiotemporal entanglement between signal and idler. 

The reduced density matrix in space and frequency for the signal, calculated by 
making a partial trace over the idler photon, writes 



P signal = Tr id i er (p) (11) 



c/q s dVt s dc\i dVtidc{ s dVt' s 
$(q s , Q s , q i; a)$*(q' s , q*, ^)|q s , f2 fl ) (q«, ^ 



and its purity is given by 



Tr (p 2 signai) = J dq,dfi s dqidfiidqidfi' s dqjdfij (12) 
x $(q s , Q s , q i; fii)$*(q' s , tf s , q*, fti) 

x$( q ;,f2;, q :,^)^(q s ,a,q^:). 



Recalling the exponential character of the mode function described by equation we 
can write equation (fT2l) as 



Tr( Psignal ) - (13) 
where C is a positive-definite real 12 x 12 matrix given by 



iV 4 exp 



-X l CX 
2 



^(q^n^qi,^ 



(14) 



x <r ( q ' s , n' s , q ,, a)$( q ;, q' s , q :, fi:)$*(q„ n s> q ;, n;). 



Notice that the difference between the matrix C and matrix B, in the last section, is 
the order of the primed and unprimed variables in the arguments of the mode function. 



Figure H] shows the space-frequency purity of the signal photon as a function of the 
spatial filter width for different values of the frequency filter width. Maximal separability 
for the photons is achieved when infinitely narrow filters in space and frequency are 
used. In the region of small values of w s and Wj, the correlations between frequency and 
space are considerable even in the case of infinitely narrow frequency filters. Different 
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Figure 4. Purity of the spatial frequency state of the signal photon for different values 
of the frequency filter bandwidth as a function of the spatial collection mode width 
(w s = Wi). The pump beam waist is w p = 400/im and the emission angle is ip = 10°. 



values for the signal photon purity can be achieved by changing the filters width, as 
it can be seen from figure HI It is important to take into account that Tr(p 2 signal ) is 
confined between the values obtained for AA S = AAj — > nm and AA S , AAj — > oo. This 
limits can be tailored by modifying the SPDC configuration. Furthermore, in specific 
configurations where non infinite filter are used in space and/or frequency, the purity 
can be tailored by using other parameters. 

Figure [5] shows Tr(p 2 signaj ) as a function of the pump beam waist w p for various 
values of the emission angle tp = ip s — In figure [5] (a) we considered a frequency 
filter AA S = AAj = 10 nm with infinitely narrow spatial filters. In this case, maximal 
purity can be found for each emission angle at a particular value of the pump beam 
waist. For a highly narrow band frequency filter with a w s = w, = 400 /im spatial filter, 
figure [5] (b) shows how the correlation between the photons is minimal for small pump 
beams being irrelevant the emission angle, for this particular crystal length. Finally, in 
figure [5] (c) the case of finite spatial and frequency filters is despicted. As can be seen, 
the purity of the signal photon is always smaller than one, it increases for noncollinear 
angles and it has a maximum for a given value of the pump beam waist. 

From the purity, equation (fl3l) . the amount of entanglement between the sig- 
nal and idler photons can be calculated by using, for example, the I-concurrence 
C = ^2(1 -Tr(p^ n J) [HUM] or the Schmidt number K = 1/Tr(p 2 signal ) [15]. 




Figure 5. Spatiotcmporal purity of the signal photon for different values of the 
emission angle, as a function of pump beam width (w p ). In (a) AX S = AA^ = 10 
nm and w s ,w, — » oo. In (b) AA s ,AAi — > nm and w s = Wj = 400 /an. In (c) 
AA S = AAj = 10 nm and w s = Wj = 400 /nm. 



5. Conclusion 

An analytical expression that fully characterizes the purity of the spatial two-photon 
state was obtained as a function of the parameters involved in the SPDC process, 
equation (jHJ). The fact that the spatial two-photon state may be mixed reveals 
a correlation between space and frequency that cannot be neglected even in the 
cases of very narrow filtering. This results were compared to previous analysis 
where the correlations between space and frequency are neglected [221 E3] and with 
other one in witch the correlations appears [5]. This correlation is mathematically 
analogous to the entanglement in the sense that the mode function $(q s , Q s , Qi, ^i) 7^ 
$ q (q s , qj)$n(^sj ^i)- However, the concept of entanglement is strongly associated to 
the possibility of a spatial separation of the correlated subsystems. 

The purity of the signal photon in space and frequency has been also calculated, 
equation (fl3j) . This expression represents a tool to tailor the space and frequency purity 
of the signal photon by means of the geometrical SPDC configuration and the filtering 
process. 
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